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Abstract 



We consider an electron coupled to a random magnetic field with local correlations and 
eventually non-zero average value. Starting from the source-term formalism of Levine, 
Libby and Pruisken, we define a generating function for static transport coefficients. We 
use a non-hermitian representation of the Hamiltonian, and a functionnal representation of 
the generating function, where the integration variable consists of non-local supermatrices, 
due to the non-locality of the Aharanov-Bohm interaction. The conductivity coefficients 
are computed using a specific saddle-point procedure, which incorporates ladder diagrams 
contributions. The situation with an external uniform magnetic field is considered with 
the two different spin projection possibilities. The limit of vanishingly small external 
magnetic field is unstable with respect to spin flip, because of E = states. In particular 
a quantized step is obtained when comparing the Hall conductance between the two 
possible spin projections. 
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1. Introduction 



The problem of an electron moving in a random magnetic field has been in the recent years 
subject of a great interest in particular because of its relevance to the Hall problem, especially 
at half filling Q. It may be defined in many ways, either on a Lattice with random fluxes 
pi |^] or in the continuum with various choices for the statistical correlations of the magnetic 
field, by considering for example random distribution of point-like vortices ]8|, |9|. The spectrum 
properties of these models are now well established. Although weak-localization has been studied 
[II], diffusion modes have been identified, and an effective field theory for weak disorder has 
been constructed fl2[| , the question of localization is however still under debate. Most recent 
results seem to converge to the idea that all states are localized except at E = for lattice 
models [13, 19]. The question of transport coefficient has been less studied. As noticed by 



Lusakowski and Turski ||20|| , there might be some additionnal features due to the random nature 
of the velocity operator. Concerning this aspect, some qualitative and quantitative perturbative 
results have been obtained in the situation of random magnetic impurities and strong uniform 
magnetic field, in particular oscillations of the Hall conductance 22, 23]. 

In the present article, we are interested in using some field theoretical methods, originally 
devised in the context of localization, to consider this question of transport coefficients. Es- 
pecially we would like, as a first step, to compute the SCBA approximation in a well defined 
manner, using this field theoretical language. We will consider a two dimensional model were 
the magnetic fluctuations are distributed with uncorrelated gaussian distribution and are added 
to an eventual uniform magnetic field. In the first part, as a preliminary, we use the source 
term formalism of Levine, Libby and Pruisken to define a generating function for static trans- 
port coefficients, which can be seen as a vacuum to vacuum transition amplitude at the Fermi 
energy, and necessitate to consider sources in matrix form. As a byproduct it is a generating 
function for Ward identities. In the second part, the model under consideration is specified, and 
a non-unitary transformation is performed, which leaves the generating function invariant and 
simplifies the averaging procedure. The generating function is then expressed, using functionnal 
integration over non-local replicated supermatrices, in consequence of the non-locality of the 
Aharonov-Bohm interaction. This is the main technical difference with the standard random 
scalar potential. The saddle-point is then exactly solved and leads to semi-circle type expression 
for the DOS. For the issue of transport coefficients, the saddle point equation is then recon- 
sidered in presence of the sources, which define in an unambiguous manner the saddle-point 
expression for these transport coefficient, and allow to eventually incorporate ladder diagrams 
contributions, which are known to be exclusively responsible of IR divergencies in absence of 
time reversal symmetry. Some perspective and natural extensions of this work, will then be 
drawn in the conclusive section. 



2. Generating function for static transport coefficients 

2. 1 Smrcka-Streda formula 

We start with the expressions for the static conductivity coefficients, obtained from the linear 
response theory. Given the general form of the Hamiltonian 

H= 2^ (p " eA)2 + V{r) (1) 
and assuming a zero temperature Fermi-Dirac equilibrium distribution, longitudinal and trans- 
verse transport coefficients can be expressed as Fermi-energy quantities, in the following way 
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a xx = —Tr[v x 8(E F - H)v x 5(E F - H)} (2) 



fxj/ — cr xy + o"^ (3) 
< = ^T^G* - H ) ~ v x 5(E F - H)G-v v ] (4) 



e 2 



< y = w Tr[(xv y -yv x )5(E F -H)} (5) 

= (E F — H db i??)" 1 are the advanced (resp. retarded) Green function. Vi = pi — eAi are the 
components of the velocity operator. V is the volume, and these expressions represent an average 
over the sample of local current-current correlation function, concerning a xx and a xy , although 
& X y is the magnetization at the Fermi energy per unit volume of the sample. It is already to be 
noticed that this last formula is not well defined unless edges are taken into account |2^, ^7|. A 
way to avoid a cumbersome regularization is to use the Streda formula pa] 



„ _ 8N(E) 



E=Ef 



(6) 



where N(E) is the integrated density of states. Since the model we have in mind is two di- 
mensional, it will be convenient to use complexe notations (z = x + iy, A z = A x + iA y , 
d z = \{d x — id y ), v z = \{v x — iv y )). Using these notations, we can recast the coefficients 
in the following manner 

2-7TP 2 

o xx = —Tr[v z S(E F - H)v- Z 6{E F - H)] (7) 
°L = % -\ 7 Tt[v z G-v- z G + -v z G + v- z G-) (8) 



xy ttV 

■ 2 

o£ = ^Tr[(zv z - zv- z )5(Ef ~ H)] (9) 



xy 

In the expression of crj , has been used the relation 



TrvzGPvzG* = --TrG ± 
2 



(10) 



a direct consequence of the commutation rules 

>,G-i] (11) 



l^G- 1 ] (12) 



2.2 Generating function 

The Field theoretic approach to localization led to realize pS| ] the intimate relation between 
d.c. transport (at T = 0) and the symmetry between advanced and retarded Green function; the 
localization regime corresponds to the restauration by quantum interferences of this symmetry 
[30, 31 1, otherwise spontaneously broken by a non- vanishing density of states at the Fermi energy. 
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An explicit illustration of this, is the source-term formalism of Levine, Libby and Pruisken [28|, 
which expresses the d.c. transport coefficients as the response of the system to some rotation 
between the advanced and retarded wave. By convenience we will reformulate slightly this 
formalism in the following. 

Consider the generating function 

J- [A z ,A g \- log _ H+ [16) 

It is left implicit here, that the operators involved in this expression acts on two-components 
wave functions, corresponding to either the retarded or advanced sector, as emphasized by 
the presence of the infinitesimal prescription irjrs. The A" are the components of the sources 
A z = A z r a , which are of matrix form, r a beeing the usual Pauli matrices 



In addition, we have 



H(A z ,Az) = Uv z - ^A z ),(v- Z - + V(z,z) (15) 



(where {,} means anticommutator). If we restrict now the sources to be spatially uniform for 
a = 1, 2 and 

iB 3 . . 

A = —z (16) 

B 5 being a uniform magnetic field then, the static transport coefficients can be expressed in the 
following way 



1 d 2 T 



2ttV dA\dA 

i - -_L d2jr 

° xy ~ ttV dAldAj 



A=0 



A=0 



xy 2inV dB 3 

by direct comparison with (0,8,9), using the expansion 



A=0 



(17) 
(18) 
(19) 



T\A a z ,A%\ = logdetfl - AHG] = -Tr[AHG + ^{AHG) 2 + ...] (20) 

/ Q+ Q \ 

with AH = H(A Z ,A Z ) — H and G = f ^_ J , and using the Ward identity (|To|) which in 

fact is a consequence of global U(l) x U(l) invariance of T [28]. The relation concerning a l J y is 
reminiscent of the Streda formula (|6|). The connection is simply made by realizing that 

dT[A%Al] 



dE F 



Tr[G + {B + B 6 ) +G~(B- B A ) - G + {B) - G~ (B)] (21) 



(where the argument of the G"s denotes the sum of the magnetic field B and the fictious one 
coming from the source A 3 ). Since 

G ± (B) = Re(B) ^ mp(E F , B) (22) 
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(Re(B) is the real part of the Green function, and p(Ep, B) is the density of states at the Fermi 
energy and for a given magnetic field B). Taking the first order term in the last relation, and 
integrating over energy leads to the relation 

d& = - 2mV —dB- (23) 

which makes the connection with the Streda formula. The underlying reasons for these relations 
concerning o~j.y is that the source A 3 is the conjugate field of the persistent current, locally given 
(at position z = (z, z)) by 

j° z (z) = ^-Tr[v z \z><z\Gr 3 ] (24) 
This entails, that if A 3 corresponds to a magnetic field, 

B 3 = l -{d z Al - d- z Al) (25) 
then this magnetic field is conjugate to the magnetization (at the Fermi-energy), locally defined 

by 

*M = sF/ A (i^>-rb**>) (26) 

meaning that this quantity can be obtain from the generating function in the folowing manner 

M « = -sbl^L (27) 



This is quite similar also to formulaes derived in the context of scattering theory |32|, 32 , 34 
where the logarithm of the determinant of the scattering matrix is used as a generating function 
instead of T. 



3. The model 

3.1 Pauli equation 



After this preamble concerning the transport coefficients, we are now in position to define 
the model we are interested in, namely to specify equation (|l|). As already mentioned in the 
introduction, the random magnetic problem should be relevant for studying transport in a 
fractionnal quantum-Hall state, especially at even denominator filling factors. In this case, the 
ground-state consists of an assembly of composite Fermions, forming a weekly interacting Fermi- 
sea, and scattered by randomly distributed impurities. Due to the Chern-Simons gauge field 
needed to describe the ground-state, the effect of impurities is to produce fluctuations in the 
fictious magnetic field by modifying the local density of particules yj . This entails the presence 
of magnetic vortices, located on the top of these impurities. For example, in a v = - fractionnal 
state, a delta impurity is expected to be accompagned with q quantum of flux. As a consequence, 
a relevant model for this problem could actually be defined by 

H=^-(p-eA) 2 +XB(r) (28) 

where B(r) is a random magnetic field plus the eventual effective external field, A the cor- 
responding static gauge field, and A a proportionality constant between the impurity scalar 
potential an the corresponding induced magnetic field. Once the statistical properties of the 
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magnetic field are specified, the model is entirely defined. Instead of studying this model for 
an arbitrary value of the parameter A, we will restrict ourselves to the two particular values 



A = ±^-. In the context of the Aharonov-Bohm problem [35, 36], or Chern-Simons theory 
these particular values correspond to the situation where the scale-invariance of the theory is 
not broken quantum mechanically and can be viewed as fixed point of the parameter A under 
renormalization. They are moreover interpreted to be the spin-coupling constant of the electron 
with the magnetic field []38|| , the equation (|28|) beeing then precisely the Pauli equation, with 
the sign depending on the projection of the spin along the z axis. In the present context, this 
particular choice has as well the virtue of eliminating ultraviolet divergences, which otherwise 
occur at arbitrary A in the perturbation expansion (of the average density of states for example), 
indicating also the scale invariance of the theory for these particular values. As a consequence 
we will consider this to be the most relevant situation, and postpone the analyses for arbitrary 
A to further investigations. 

For simplicity, we will consider a locally correlated random magnetic field, in its most simple 
expression, namely the white noise case, 



B(t)B(t')=w6(t-t') (29) 

and with all the higher order cumulants set to zero. If we compare to a poissonian distribution 
of 5 fluxes (f), with mean density n, the white noise distribution amounts to conserve only the 
second cumulant with a coefficient equal to 

w = ncj) 2 (30) 

This approximation, to some extent, is justified from the field theoretical point of view, for 
scalar impurities (not magnetic). The second cumulant gives rise to a local quartic interaction 
term, and other cumulants contribute to terms of higher degre in the interaction, which are 
in principle irrelevant. This might be less justified in the random magnetic-field problem (see 
reference || for the average density of states for poissonian magnetic impurities, M for non- 



locally correlated poissonian magnetic impurities and [22] concerning the average conductivity 
for poissonian magnetic impurities), where as we shall see, the interacting term hapens to be 
non-local. 



3.2 Non-unitary transformation 

The model beeing now defined let us consider it more explicitely. Given a magnetic field 
distribution B(z), the static gauge field can be expressed in the Coulomb gauge (divA = 0), 
using complex coordinates, according to 

(recall that z means (z, z), B(z) is not supposed to be holomorphic) . In the following we will 
have to distinguish between two situations, whether, in addition to the fluctuations, there is or 
not a uniform magnetic field. 

When the average value of the magnetic field is non zero, then the two possible signs for 
the spin projection are unequivalent, and assuming by convention this uniform field to have the 
positive value Bq, along the z axis, the Hamiltonian takes then the two possible forms 

H± = -2 (d z - JL J d V ^S(«0) [d- z + i- J dV j^Btf) + (1 ± l)S(z') (32) 
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Where the value of e and m are from now on set to 1. Having in mind to average later on over 
disorder, we perform a non-unitary transformation on the wave-function, in order to eliminate 
the quadratic term in the fluctuating magnetic field || 

V,( z )= e ±2¥ J" d 2 z6(z) log |*-*'| ^( z ) (33) 

b = (B — Bq)/Bq is the dimensionless fluctuating magnetic field. The generating function T 
which will be used later, is invariant under this transformation, as long as the Hilbert space 
is not modified. This is indeed the case, qualitatively because 6(z) average itself to zero over 
the plane, consequentely wave functions remain square integrable under this transformation, or 
under the inverse. At this point it is convenient to introduce the raising and lowering Landau 
level operators, 

a = d- z + ^z (34) 

a f = ~d* + \z (35) 

(the cyclotron energy is set to one) The Hamiltonian acting on the wave function ■0 takes 
then the simple form 

H + = 2a«t + - [ dV-^V (36) 

7T J Z — Z' 

H- = 2a)a-- [ dV ^-a (37) 



This last expression, corresponding to the spin down case, makes the (expected) zero mode 
directly apparent. It consists of the lowest landau itself, giving therefore a degeneracy to these 
particular states in accordance with the Atiyah Singer theorem. The fact that this feature 
explicitely appears here, is due to the very analogous form of the factor multiplying tp with the 
Aharonov-Casher solution [39]. 



In absence of a uniform magnetic field, considering the statistical properties of the random 
field, there is no priviledged orientation given to the plane. The sign of the spin projection 
is therefore indifferent, and taking it negative by convention, we then perform the same non- 
unitary transformation (minus sign in the exponent). The Hamiltonian acting on ip has then 
the expression 

H = -2d z 3- z + - [ d 2 z'^-d z (38) 
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4. Saddle point DOS 

4-1 Replica Super symmetry representation 

We turn now to the functionnal approach to the question of computing the d.c. transport 
coefficients. The objective is to average the generating function over disorder, which necessitates 
to express it in a convenient form. Although we know that im principle, in order to compute the 
mean conductivity, it is sufficient to average exp(J r ) (the equivalent of the partition function, 
instead of the free energy), we will follow the logic of this presentation, and average T itself. 



Using the combination of the supersymmetry method of Efetov [31] and the replica trick, it is 
possible to express T in a suitable form for average. 
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The supersymmetric calculus gives the possibility to view a ratio of two determinants as a 
superdeterminant 

det[E F -H(A z ,A s ) + ir]T3] f . ... , f . 

— — — — = = Sdet[E - H{A Z ,A Z ) + ir]T 3 \ (39) 

det[Ep — H + irjT 3 \ 

where the numerator corresponds to the grassmann sector and the denominator to the bozonic 
one. The sources are still ordinary complex fields, but are acting in the Grassmann sector. The 
superdeterminant is then expressed in terms of a Gaussian integral over superfields. Moreover, 
the replica-trick amounts to express the logarithm according to 

logSdet[E-H(A z ,A 2 )+irjT 3 ] = lim —Sdet[E — H(A Z ,A Z ) + irjT 3 ] q (40) 

<?^o dq 

which has to be understood by means of analytical continuation of q to real values. This allows 
to express J- in the way 

F[A?,A%] = lim 4- / DVD* e iJ^HE-H(A z ,A,)+irir 3 xl q ^ (41) 
q^o dq J 

where the integration is performed over 2q components superfields, q complexes and q grassman 
components. In this form it is now possible to average over the magnetic fluctuations. Let us 
consider the case with no uniform magnetic field. The term to be averaged is 

< e -^/d 2 zdVf(z)^i3^(z) >b = e ^/dWz2*(z 1 ).D* 1 *(zi)§^g *(z 2 )D Z2 *( Z2 ) (42) 

where D z replaces d z because of the sources and where has been utilized the identity 

d 2 z' _ 1 = (43) 

(zi - z')(z 2 - z') Z1-Z2 

As already anticipated, the field theory we are obtaining, contains, in addition to the free part 
of the action, 

So = i J d 2 z*(z) (E + 2D Z D- Z + i V T 3 x /,) (z) (44) 

a quartic interacting term, which in contrary to the usual random scalar potential, happens to 
be non local, because of the long range nature of the Aharonov-Bohm interaction. This term 
can be reexpressed using the definition of the supertrace 



w 

Si n t = — I d zid z 2 Str 



^—^D Zl *{z l )*{z 2 )D Z2 *{z 2 )*{z l ) 

z\ — z 2 



(45) 



and can nevertheless be decoupled, using a Hubbard-Strattanovich transformation, that is to 
say the identity 



c f — -^-StrQQ +[ d 2 z 1 d 2 z 2 Str 



£)z 1 *(zi)*(z2)Q(z2 1 zi)+Q(zi,Z2)D Z2 *(z 2 )*(z 1 ) 



e *im = J DQDQ e w J ^ l ~" 2 1 ' J (46) 

The functionnal integration involve 2nx2n non-local supermatrices with no particular constraint, 
in contrast to the local theory |HJ. Permutating again the supertrace in Si n t, and integrating 
over the 2n components superfields, we obtain 



d f — Str log 



E+2D z D E +irjT 3 xI q +i^-^QD z +iQD 



^StrQQ 



F[A°,A?] = limj- j DQDQ e L n-*,- - -j » " ( 47 ) 
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4-2 Saddle point approximation 



The saddle point is then obtained as the stationnary point when Q and Q are varied inde- 
pendantely. Expanding at first order the supertrace of the logarithm in the above expression, 
leads to the set of equations 

w 

<3o(zi,z 2 ) = i— < zi| d z G \z 2 > (48) 

Z7T 

Q ( Zl ,z 2 ) = i^- < zi| d z G |z 2 > (49) 

Z7T Z\ — Z2 

G (zi,z 2 ) = < zi| — 1 |z 2 > (50) 

£ + 2d z d g + 2tQ d z 

the sources beeing set to zero for the moment. The main difference with the scalar case is that 
these relations are equalities between operators. In absence of a uniform magnetic field, we are 
seeking a solution which preserves the translationnal and rotationnal invariance, meaning that 
G has to be a function of p = |p|. In addition we look for this solution in a diagonal form. In 
that case the equation for any diagonal components g(E,p) of G(p) rewrites 

9( E ,P) = p — l w 2 (tp \ ( 51 ) 

This equation admits two real solutions 

9 ± {E,p) = ^-Ae-\p>± ^{E-Ev){E-E+)) (52) 
wp \ z / 

outside the energy interval [E~,E+] with 



One corresponds to a stable saddle point, the other to an unstable one. By looking at the 
sequence g n+ \(E,p) = l/(E— \p 2 — j^p 2 g n (E,p)) one observe that when E < E~, the attractive 
fixe point (stable) is g + (E,p), even though it is g~(E,p) when E > Ep . This is then sufficient 
to analytically continue g(E,p) onto the complexe plane, without any ambiguity, 

2vr / 1 



9( Z 'P) =^{ Z ~2 p2 ~ V {Z ~ E ^ {Z " E ^J (54) 

where the square root is defined by using the convention that the argument of Z — E~ and 
Z — E+ vary in the interval [— n, tt]. As a result we have a branch cut situated on the segment 
Z€[E~,E+], 

g(E ± in,p) = ^(e-\ P 2 t iyj (E+ - E)(E - Ef)) (55) 
wp \ z / 

corresponding to a non-vanishing contribution to the density of states 



l>,m - ^: (Gp(E ~ in) - G P (E + in)) = ^ \/( E p ~ E ){E - E p ) (50) 



Therefore, the SCBA density of states is obtained by summing this last expression with the 
measure pdp/2ir, which leads to 
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p(E) = \±(2E + f) ~%<E<0 (57) 

lap 

and is responsible for the spontaneous breaking of the symmetry V(2n) — ► F(n) x F(ra), i.e. 
between the advanced and retarded sector, at the saddle-point 



2vr /„ 1 2 



G (£,p) = — [E - -v l - ir 3 x / gV /(i?+ - £)(£ - £ p ~) (58) 

The negative part of the DOS is actually unphysical, because the spectrum of this problem is 
supposed to be positive. The possible origin of this indesirable feature shall be discussed in the 
conclusive section. 

4-3 Broadened Landau spectrum 

If we turn now to the case, where a uniform magnetic field is present, we are expecting to 
find a saddle-point solution which is both invariant under magnetic translation and rotations 
around the origin (if we choose the symmetric gauge centered at the origin) . The representation 
of this group of transformation is irreducible in each Landau level, meaning that G(E) has to 
be decomposed on the set of projection operator P n corresponding to each Landau Level. In 
other words, G(E) has to be a function of aa) 

oo 

G(E,aa^) = J2 G n(E)P n (59) 

n=l 

For the spin-up case, the saddle point equation (with the sources set to zero again) takes the 
form 

G(z,,z;) =< zi| ^ _ _ I j-^ , G t N > («0) 

7T Z — Z 

(where a shorthand notation has been to express that the kernel of G a) is multiplied by 
jEjr)- Using the properties (see appendix) that 



zi - z 2 



< zila 1 " G(£',aa t ) a t |z 2 >=< zi|aa t G(S, aa t )|z 2 > (61) 



and following the lines of reasoning described in the preceeding case, we obtain (with n labelling 
the set of Landau levels) 

G (E, n) = it 2 (e n (E) - ir 3 x I qPn {E)) = = (62) 

E - 2n + ir 3 x I q \J {E - En)(E+ - E) 

with 

E± = 2n± Z^^n (63) 
which leads to the standard semi-circle law for the SCBA density of states 



^/Aan - (E/b - 2n)' 2 E~ < E < E+ 

otherwise 
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where the parameter a = ^- has been introduced, and where the dependancy with the b field 
has been restored. In this form it is explicitely seen that p(E), and actually any other quantity 
is a function of E/b and a, the scaling property of the system is not broken by any UV cut-off. 
The real part of the Green function is given by 



e n (E) 



1 



2ir 2 na 



' E/b -2n + yJ{E/b - 2n)' 2 - Ana 

E/b-2n 
t E/b -In - yJ(E - 2n)' 2 - Ana 



E < E~ 



E- <E<E+ 



(65) 



E>Et. 



The DOS is therefore composed of a sequence of enlarged Landau levels, with a width increasing 
in proportion to the square root of the LL index n. The result is not very different from what 
is obtained for a poissonian magnetic impurities system |8|, This is however different from 
the scalar random potential situation, where the width is not n dependant, and this fact has a 
simple semi-classic interpretation f4C| , ||; the semi-classic orbits in the Landau problem consists 
of cyclotron orbits, the number of time the electron circulates around such an orbite corresponds 
to which Landau level it belongs to. Moreover, the magnetic flux passing through a given area, 




Figure 1: Broadening of the Landau Levels at fixed b for the values a = 10 2 ,310 2 ,10 1 , 1, 10 of the 
disorder parameter. 

has fluctuations proportionnal to the square root of this area, meaning that the fluctuations of 
the flux passing through a cyclotron orbit surrounded n times, is proportionnal to y/n, giving the 
typical width of the broadened Landau levels. In the spin-down case, the saddle point equation 
reads 



G(zi,z 2 ) =< zi 



E-2a\a-^~^ 1 a G a 



Z2 > 



and using now the property 



zi - z 2 
z\ - z 2 



< z\\a G(E,aa^) a\z2 >=< zi\a^ aG(E, aa' )\z 2 > 



(66) 



(67) 
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we obtain the same expressions but with the modification n 
Landau level at E = remains unperturbed. 



n — 1. In particular, the lowest 



5. Saddle point Conductivity 

5. 1 In absence of a uniform magnetic field 



The next step is to calculate the transport coefficient at the saddle point. Using the gen- 
erating function, it amounts to replace the integration over the Q and Q matrices in (^) by 
the saddle-point solution Qq and Q Q of equation fl48|). Since the saddle point doesn't break the 
replica symmetry, we therefore obtain the saddle point expression for the generating function 



Tq[A^A%\ 



logSdet E + 2D z D z + 2iQ D z 



(68) 



using the fact that | L E^ Qo( z i, z 2) = Qo( z i' z 2)- By differentiation with respect to the sources, 
according to receipes given in part one, we would then obtain bare values for transport coef- 
ficients. However, at this point, we have to realize that the saddle point generating function 
doesn't fulfil local U(l) x £7(1) (gauge) invariance, namely because Qq and Q are non-local op- 
erators. This originates in their dependancy with respect to the operators d z and d z accordingly 



to equation (48). A straightforward procedure to restore gauge invariance would be then, to re- 
place these operators by D z and D z . However, this procedure is ambiguous because T3 does not 
commute with D z and D z , leaving any expression of Qq to be arbitrary. Instead, if we compute 
the saddle point, while keeping the sources fixed, we obtain the correct dependancy of Qq and 
Qq with the sources. Let us denote Q(A) and Q(A) the corresponding saddle point solutions. 
From the first part we know that it is sufficient to consider uniform sources to compute the 
conductivity tensor. In that case the saddle point expression for the generating function can be 
written in the form 



Fo[A%,Af] = logSdet [E + {D Z ,D- Z } + iQ(A)D z - W- Z Q{A) 
with the saddle point constraint 

Q(A) 
Q{A) 

G(A) 



—StrQ(A)Q{A) 

w 



i^-D z G(A) 
G(A)D~ Z 



E + {D z ,D- z } + £ {G(A),D- Z D Z } 



(69) 

(70) 
(71) 
(72) 



The computation of transport coefficients necessitates then to solve this last equation. Assuming 
that G(A) is analytic with respect to the sources, it is then sufficient to solve at first order in 
the source. Consider the expansion 



G{A) =Gq + G z a A a z + GlA% + .. 



(73) 



Deriving (p9|) twice with respect to the sources, according to the saddle point constraints leads 
to the following expression for the conductivity coefficients 



&XX 

I 



1 + 



a 



■>-y 



AttV 

i 

2fnV 



w , + 

Vz + —{VzGq + GqV, 



47T 



■w 



(74) 
(75) 
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where the trace is yet restricted to the Grassmann sector, and Gq on this sector reads 



whereas 







o Gi 



G + 






G n 



(76) 



(77) 



In absence of a net magnetic field, we expect G z a to be a function of p and p because of transla- 
tionnal invariance. Solving ([7^) with this assumption we obtain 



P 
~1 



1 + M.(G + Gl)) Go Gl 



1 



^P 2 G Gl 



a 



1,2 



(78) 



i- 



the striking feature of this last expression is that it is not defined when p is such that E p < 
Ep < the denominator vanishes because in this interval we have (fjq) 



GqGq 



47T 

wp 2 



(79) 



The reason for this is that the procedure of keeping sources in the computation of the saddle 
point, amounts to incorporate ladder diagrams contribution, when the trace is taken. These 
diagrams are in principle responsible for weak localization corrections, and since the system, 
because of the magnetic fluctuations, is not time-reversal invariant, they are actually the only 
ones. Therefore this divergency indicates that G(A) is not an analytic fonction with respect 
to the sources, which has to be related to the eventuality that the system is localized. The 



expression of G*(p) suggests the possibility of a self-consistent treatment of the problem [41|; 
this would however require a more complete treatement, beyond saddle point approximation, 
which is out of the scope of this paper. 

5.2 Conductivity in presence of a uniform magnetic field 



The situation is quite different when an external uniform magnetic field is added to the 
problem. In that case, the IR divergency observed preceedingly happens to be removed. Indeed, 
the solution to ( |72| ) can be searched in the form 



G z a (a,a)) = a/g(aa>)T a 

and using the rules (for an arbitrary function /), 

af{a [ a) = f{aa))a 
oVCoot) = f{a)a)a ] 

we obtain the solution 



G*(o,a t ) = -i 
for a = 1, 2 and 

Gf(a,at) = -i 



1 + ^ (Go(aat) + Gj(ata))] G Q {aa)) G^a) 



1 - ^(2aat - 1) Go(aat) Gj(ata) 
1 + £ (Go(oot) + G (o+o))] G (aat) G (a!a) 



a' r n 



(80) 



(81) 
(82) 



(83) 



1 - ^(2aat - 1) G (oat) G (at a ) 



a 1 " t 3 



(84) 
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Figure 2: Transverse conductivity in unit of %- for the spin-up (left) and spin-down (right) cases for 

a = Kr 2 ,icr' 

filling number. 



a = 10 2 , 10 1 , 1, 10, compared to the classical value (straight line). The x axes is the Landau level 



where Go(aa^) is defined according to ( ]59| ) and (|62]). Inserting this last expression in ( p7| ) and 
(|89| ) allows to compute the conductivity tensor. However for a xx it is convenient to use the 
following Ward identity which is a consequence of global U(l) x J7(l) invariance, 



Vz + — {Vz&o + G V 2 



Gfr 3 } = (85) 



and which is the saddle point version of fllQ|). Defining 

g s n = ir 2 (e n (E F ) - is Pn (E F )) (86) 



with s = ± and e n (Ep) and p n {Ep) given by (p5|) and ( p4[ ) we can then express a xx in the form 
(recall a 



w > 



= E E ss n - 1 ^-j — ^ — t^~tt> — ( 87 ) 

In the strong magnetic field limit, i.e. when the Landau Levels do not overlap, it is quite easy 
to see that o xx vanishes. Surprisingly enough this hapens to be true whatever b is (numerical 
check) 

&xx = 

Concerning the transverse conductivity, we obtain for the dissipative part 



n=l s=± 2 



although the contribution comming from the magnetization is obtain by using the Streda formula 
© 

1 00 1 

41 = l N( - E ^ ~ E 2 {Ep/b + 2n )Pn( E F) (90) 

n=l 
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Altogether the sum simplifies exactly to 



up 
'xy 



n=l 



p(Ef) 



(91) 



For the spin-down case, all the n appearing in the fraction of expression ( |87| ) and ( p9| ) are shifted 
by —1. The conclusion concerning a xx is not modified, the longitudinal conductivity vanishes 
exactly for all values of the magnetic field. For a xy we obtain a slightly modified expression, 
exept for the sign relative to the "classical" value, 



a 



„down 

xy 



-N(Ep 



OO 771 
I \ " ^ F 



(92) 



n=l 



which again is obtained with no approximation. The Hall conductance is plotted in figure 2 
for the two orientations of the spin. In the strong magnetic field limit we obtain oscillations of 
the transverse conductivity, in qualitative agreement with the conclusion of reference |2^] for 
the spin-down case. The spin-up plot ressembles qualitatively the numerical results obtained in 



reference [42]. 




0.5- 



'-V' "-!>" '"V " 'i " "4' " '& " "4" 1 'ib' ' ' 'V ' 

E 

2 

Figure 3: Difference between the spin- up and the spin-down transverse conductivity in unit of %- for 
b = 5, 1, .1, .01 and w = ir, versus Fermi energy. 

The limit of vanishing mean magnetic field is a bit more mysterious. It should be first noticed 
that function of electronic integrated density is divergent. Although paradoxical, one 

can consider this to be coherent with the fact that a xx vanishes, if one think of the semi-classical 
relation between these two coefficients ]|43|]. Anyway, something interesting happens when we 
compare the Hall conductance for the two different spin projections in this limit, as a function 
of the Fermi energy (see figure 3). This difference takes the quantized value of 2 (in unit of 

2 

%-) when b goes to zero. The first quantum is simply traced back to the additionnal Landau 
level lying at E = in the spin-down situation and which is not affected by disorder (see the 
strong magnetic field curve b = 5 in figure 3). The second one seems to indicate that the system 
restores the chiral symetry (considering spin-up and spin down situations as a whole) when 
the average b goes to zero, by reintroducing progressively a zero energy state in the spin-up 
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spectrum. In addition this state restores parity, by contributing with the opposite sign to the 
Hall conductance. In this sense, the limit of zero magnetic field is unstable with respect to 
spin-flip at E = 0. 



6. Conclusion 

The formalism presented in this paper, which proposes to compute directly static transport 
coefficients from a generating function, has been used at the basic level of the saddle point 
approximation. Already at this stage was exemplified the possibility of incorporating first dom- 
inant contributions, namely the ladder diagrams contribution, by keeping uniform sources when 
computing the saddle point. In the exemple we considered, the white noise random magnetic 
field, it was observed that in absence of an external uniform magnetic field, the saddle point 
solution is not analytic with respect to the sources indicating a non-simply diffusive behavior. 
Slowly spatially varying sources should be considered in order to analyse further this question, 
and a self-consistent argument analogous to the one of Vollhardt and Wolfle for the scalar disor- 



der case 1 41] could be elaborated. Furthermore, this formalism should be possibly generalized for 
the computation of dynamical transport coefficients, and eventually extended to the non- linear 
regime. 

Concerning the random magnetic field problem, the free parameter A which was introduced 
in section 3, was only considered for two particular values, which presumably insures the system 
to remain scale invariant quantum mechanically. The fact that the longitudinal conductivity 
was suprizingly found to remain exactly zero, when Landau levels do overlap, although not 
really conclusive (within the saddle-point approximation), is in favour of this hypothesis. Inter- 
estingly enough, although may be more technically involved, would be to consider this model 
for arbitrary A, in particular to see wether it is possible to find at this level of approximation 
a relation between a xx and a xy analogous to the one obtained by Ando et al. for the scalar 
random potential ||4q[ . In addition a flow between these two points under scale transformation 
might be expected. Beforehand the problem of undesirable negative energy states which ap- 
pear in zero magnetic field has to be cured, may be by avoiding the non-unitary transformation, 
which presumably causes some non-normalizable states to re-enter the Hilbert space. This could 
possibly clarify the nature of the zero energy state in the averge magnetic field limit, which, 
as suggested by the analyses of figure 3 seems to spontaneously break parity and to be able to 
carry one quantum of Hall current. 
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Appendix 

In this section we justify the use of the formula 

^ < zi| a* /(aa t )|z 2 >=< zi| f(aa^) a |z 2 > (93) 

z\ - z 2 
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or equivalentely 

^ -3- < zi| f(aal) a) |z 2 >=< zi| a/(aa t ) |z 2 > (94) 

where / is an arbitrary function of aa) , which is actually defined by the values f n it takes on 
each Landau Level 

oo 

/(oat) = £ / n P n (95) 

n=l 

To this end, it is convenient to introduce the raising and lowering operators of angular momentum 

b = d z + h (96) 
& f = -dz + \z (97) 

and a generating function for the set of projection operators on the Landau Levels, namely the 
thermal propagator 

oo 

G t = e~ aaH = e ' ntp n (98) 

n=l 

Then, according to the commutation rules between the a, a) and b, ¥ operators we have the 
property 

[z,G t ] = [a + b\G t ] = [a,G t } = (l-e t )aG t (99) 
as a consequence, it turns out that 

< zi| aG t |z 2 >= - _ — ( Zl - z 2 ) < zi| G t |z 2 > (100) 
In the same way we have 

< zi| a)G t |z 2 >= - ^ (zj - z 2 ) < zi| G t |z 2 > (101) 

which leads to obtain 

zi - z 2 

Zl - Z 2 

or equivalentely 

zi - z 2 



< zi| a) G t |z 2 >=< zi| G t a |z 2 > (102) 

< zi| G t a 1 " |z 2 >=< zi| a G t |z 2 > (103) 

zi - z 2 

Gt beeing a generating function for the Landau Levels projection operator, the formula is there- 
fore valid for an arbitrary function of aa) as stated above. 
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